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Abstract 

We examine various direct and indirect constraints on the lepton-specific two-Higgs doublet 
model and scrutinize the property of the Higgs bosons in the allowed parameter space. These 
constraints come from the precision electroweak data, the direct search for Higgs boson, the muon 
anomalous magnetic moment, as well as some theoretical consistency requirements. We find that 
in the allowed parameter space the CP-odd Higgs boson A is rather light (rriA < 30 GeV with 
95% possibility), which is composed dominantly by the leptonic Higgs and decays dominantly into 
t + t - ; while the SM-like Higgs boson h (responsible largely for electroweak symmetry breaking) 
decays dominantly in the mode h — > AA — * At with a large decay width, which will make the Higgs 
discovery more difficult at the LHC. Whereas, this scenario predicts a branching ratio Br{Z — > 
t + t~A) ranging from 10~ 5 to 10 -4 , which may be accessible at the the GigaZ option of the ILC. 

PACS numbers: 14.80.Cp, 12.60.Fr 
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I. INTRODUCTION 



The phenomenological success of the standard model (SM) has significantly limited the 
possibility of new physics except for the Higgs sector which remains untested. There are 
numerous speculations on the possible extensions of the Higgs sector, among which the 
simplest is to introduce one more Higgs doublet. Compared with the SM, such simple two- 
Higgs doublet models usually have much more complicated Higgs phenomenology. In the 
SM a single Higgs doublet is responsible for the electroweak symmetry breaking and the 
Higgs couplings with fermions and gauge bosons are completely determined by their masses, 
and therefore there is little guesswork in determining the discovery channels for the Higgs 
boson |]|. In the two-Higgs doublet models, however, the addition of new scalars and the 
modification of the Higgs interactions will significantly complicate the Higgs discovery at the 
LHC [2]. Given the imminent running of the LHC, the phenomenological study of various 
such models is urgently important. 

In this paper we focus on a special two-Higgs doublet model called the lepton-specific 
two-Higgs doublet model (L2HDM) {3]. Since this model is arguably well motivated from 
some fundamental theory and also has some phenomenological virtues (e.g., it can provide 
a natural explanation for the leptonic cosmic ray signals reported by PAMELA and ATIC 
4]), it has attracted much attention . We will check various constraints on the model 
parameters and then scrutinize the property of the Higgs bosons in the allowed parameter 
space. These constraints come from the precision electroweak data, the direct search for 
Higgs boson, the muon anomalous magnetic moment, as well as some theoretical consistency 
requirements. Our main observation is that in the allowed parameter space the CP-odd 
Higgs boson A must be light (rriA < 30 GeV with 95% possibility ), which is composed 
dominantly by the leptonic Higgs and decays dominantly into rr; while the SM-like Higgs 
boson h (responsible largely for electroweak symmetry breaking) decays dominantly in the 
mode h — > AA — > 4r with a decay width usually exceeding several tens of GeV, which may 
make the Higgs discovery more difficult at the LHC. 

This paper is organized as follows. In Sec. II we recapitulate the L2HDM model. In Sec. 
Ill we examine various constraints on the parameter space and study the properties of the 
Higgs bosons in the allowed parameter space. Finally, In Sec. IV we give our conclusion. 
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II. THE LEPTON-SPECIFIC TWO-HIGGS DOUBLET MODEL 



The L2HDM is a special two-Higgs doublet model in which one Higgs doublet <j>x couples 
only to leptons while the other doublet (fi 2 couples only to quarks. Both Higgs doublets 
contribute to the electroweak symmetry breaking: v 2 = v 2 + v \ = (246 GeV) 2 with v\ and 
t>2 being respectively the vacuum expectation values of <\>\ and <p2', whereas their relative 
contributions can be quite different and can be parameterized by the ratio tan/3 = i>2/Vl- 
So for a large tan (3 the lepton Yukawa couplings can be greatly enhanced. 

The Yukawa interactions and the Higgs potential are given by 

c Y = -YiHifaj - v;vv', - v;v>>r/ ; + h.c. (i) 

+X 3 \h\ 2 \h\ 2 + A4|<M02| 2 + y [(0l02) 2 + h.c] , (2) 

where i,j are generation indices, Y e , Y u and Yd are 3x3 Yukawa matrices, q{ and li denote 
respectively the the left-handed quark and lepton fields, U{ and d{ denote respectively the 
right-handed up- and down-type quark fields, denotes the right-handed lepton fields, and 
m 2 and A are free parameters. 

Just like the usual two-Higgs doublet model [3] , the spectrum of the Higgs sector includes 
three massless Goldstone modes, which become the longitudinal modes of W and Z bosons, 
and five massive physical states: two CP-even states h and H, a pseudoscalar A, and a pair 
of charged states H ± . These states are related to the doublets 4>i and (f>2 by 

0° = - ^(^i + Hcosa — hsina + iG° cos (3 — z'Asin/3), (3) 
v2 
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(v 2 + H sin a + h cos a + iG° sin (3 + iA cos 0), (4) 



0± = G ± cos/3-F ± sin/3, (5) 
0± = G ± sin /3 + H± cos /?, (6) 

where a is the mixing angle that diagonalizes the mass matrix of the CP-even Higgs fields. 

Due to the constraint v\ + v \ = (246 GeV) 2 , the eight free parameters in Eq. (j2J), i.e., 
Aj (i — 1, . . . , 5) and m 2 (i=l,2,3), reduce to seven. In our analysis we choose the following 
seven parameters as the input parameters of the L2HDM: 

rrih, rriH, ^a, %±, tan/5, sin a, A 5 , (7) 



where rrih, tua,, tuh, and are the masses of the corresponding physical states. Through- 
out this paper, we use H (h) to denote the Higgs boson with (cj)®) as its dominant 
component, which means that we choose cos 2 a > 1/2. 



The interactions of the Higgs physical states with fermions are then given by 

Ly = -(cos a e^-H — sin a e^n) 

2m w cos p 

gm Q . 

(sin a qiqiH + cos a qiqih) 



6] 



2mw sin (3 

igm u . _ ignid - 

+- cot/? ua 5 UiA - cotp dijsdiA 

2m w 2m w 

igm e . 

+- tan/3 e^eiA 

2m w 

+ J^ i3 cot/5 Ui{m u .P L - m dj P R )djH + 

+-^^ta,nPi? i P R e i H + . (8) 
y2m w 

Obviously, for a large tan j3 the lepton Yukawa couplings are greatly enhanced relative to 
the SM prediction. One can also check that the couplings of ZZh and ZZH are given by 

V ZZh = -^f-anOg-a), (9) 
cos t) w 

Vzzh = cos(/3 - a), (10) 

cos W 

which satisfy the sum rule V§ zh +Vz ZH = Vz Z h SM - ^or a large tan (3 (this is the case required 
by the experimental constraints, as shown below), the coupling ZZh is dominant over ZZH , 
so h is usually called the SM-like Higgs boson. 



III. CONSTRAINTS ON THE L2HDM 

We note that both the theoretical consistency and the electroweak data have limited 
the parameter space of the L2HDM. In our study we consider the following theoretical 
constraints: 

(1) The perturbativity is valid in the Higgs sector, which requires Aj < 47r (i — 1, . . . , 5). 

(2) The 5-matrix satisfies all relevant tree-unitarity constraints, which implies that the 
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quartic couplings Aj satisfy 

3(Ai + A 2 ) ± v / 9(A 1 -A 2 ) 2 + 4(2A 3 + A 4 |) 2 < 16tt, 



Ai + A 2 ± v/(Ai -A 2 ) 2 + 4|A 5 | 2 < 16tt, 
Ai + A 2 ± - A 2 ) 2 + 4|A 5 | 2 < 16tt, 
A 3 + 2A 4 ±3|A 5 | < 8tt, 
A 3 ± A 4 < 87r, 
A 3 ± |A 5 | < 8tt. 



(11) 



(3) The scalar potential in Eq. (jSJ) is finite at large field values and contains no flat 
directions, which translate into the bounds [7J] 

Ai, 2 > 0, 

A 3 > -\A1A2, 

A 3 + A 4 - |A 5 | > -aAiA^. (12) 

On the experimental side, we consider the following constraints: 

(4) The lower mass bound on the charged Higgs bosons: > 92 GeV 

(5) The constraints from the LEP search for neutral Higgs bosons. We compute the 



signals from the £ 



26, 2r, 46, 4r, 262r 



iggsstrahlung production e + e" 
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> (fli = /i, H) with fli -> 
llj and from the associated production e + e _ — > ii/jA with 



ifj-A — > 46, 4r, 262r, 66, 6r [12| and compare them with their LEP data. We also 
consider the constraints from e + e~ — ► ZiJj by looking for a peak of M#. recoil mass 
distribution of Z-boson [3] and the constraint of Y(Z — > < 5.8 MeV when 

m A + m Hi < m z [14]. 

(6) The constraints from the LEP search for a light Higgs boson via the Yukawa process 
e + e~ — ► ffS with f — b,r and S denoting a scalar 151 ] . These constraints can limit 
ffS coupling versus m>s and thus can constrain the parameters of the L2HDM. 

(7) The constraints from the H^-boson mass. The L2HDM Higgs sector can shift the W- 
boson mass through radiative corrections. We require the corrected H^-boson mass to 
lie within the 2a range of the global-fit value. The SM prediction for the H^-boson 



mass is 80.363 GeV for rru = 173 GeV and m# = 1 



1 GeV [16j, and its fitted value 



is 80.398 ± 0.025 GeV 17]. We use the formula in 18| in calculating the mass and 
consider the effect of a different top quark mass (in our calculation we take m t = 171.3 
GeV). We also subtract the contribution from the SM Higgs boson to avoid double 
counting the contribution from the Higgs sector. 

(8) The constraints from Zt + t~ coupling. For a large tan/3 the L2HDM Higgs sector can 
give sizable radiative corrections to Zt + t~ coupling. We calculate such corrections 
and require the corrected Zt + t~ coupling to lie within the 2a range of its fitted 
value. The SM prediction for this coupling at Z-pole is given by gy M = —0.03712 and 
gSM = _o.50127 [16J, and the fitted value given respectively by -0.0366 ± 0.00245 



and -0.50204 ±0.00064 



16]. We use the formula in 18[ in our calculation. 



(9) The constraints from r leptonic decay. We require the L2HDM correction to the 
branching ratio Br(r — > eu e u T ) to be in the range of —0.80% ~ 1.21% We use 



the formula in 



19] in our calculation. 



(10) The constraints from the muon anomalous magnetic moment a^. Now both the theo- 
retical prediction and the experimental measured value of a M have reached a remarkable 



precision, but a significant deviation still exists: a e * p — a^ M = (29 ± 8.8) x 10~ 10 |20| 
In our analysis we require the L2HDM to account for such difference at 2a level. Note 
that in the L2HDM, gets additional contributions from the one-loop diagrams in- 
duced by the Higgs bosons and also from the two-loop Barr-Zee diagrams mediated 



by A, h and H 21]. If the Higgs bosons are not too light, the contributions from 
the Barr-Zee diagrams are more important. To account for the discrepancy of a M , one 
needs a light A along with a large tan f3 to enhance the effects of the Barr-Zee diagram 
involving the r-loop. The CP-even Higgs bosons are usually preferred to be heavy 
since their contribution to is negative. 

11) Since the CP-odd Higgs A can be quite light and h, H — >• A A may open up with a large 
decay width, we require the width of any Higgs boson in the L2HDM to be smaller 
than its mass (otherwise the Higgs boson may be too fat). 
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With the above constraints, we scan the parameter space of the L2HDM in the ranges: 

l<tan/3< 200, 5 GeV < m A < 100 GeV, 

5 GeV < m KH < 350 GeV, 92 GeV < m H+ < 350 GeV, 

-y/2/2 < sin a < y/2/2, |A 5 |<4tt. (13) 

150 , , , , , , , , , 




40 



(GeV) 



FIG. 1: Scatter plots of the allowed parameter space in the plane of tan/3 versus tua- 



With 10 12 random samplings, we get the allowed parameter space shown in Figs. 1-2. Fig.l 
shows that the allowed parameter space has a light A ( uia ^ 80 GeV) and a large tan (3 ( 
37 & tan (3 ^ 130), which mainly comes from the explanation of the a M discrepancy. Among 
the surviving samples displayed in Fig.l, about 95% satisfy uia < 30 GeV and about 70% 
satisfy uia < 20 GeV, which means that a very light A is highly preferred by the constraints. 

Fig. 2 shows the allowed parameter space projected in the planes of m h and m H versus 
rriH+. Three characters should be noted about this figure. The first is that all the Higgs 
bosons are lighter than 350 GeV (lighter than 250 GeV for about 90% of the surviving 
samples), which is mainly due to the unitary requirement and the constraint. The second 
is that h and H can be as light as 58 GeV because the LEP2 bound is relaxed significantly 
due to the weakened ZZh and ZAH couplings by the seizable mixing angle a and the open- 
up of the new decay mode H,h — > AA — > 4r 29]. The third character is that the values of 
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FIG. 2: Same as Fig.l, but projected in the planes of and mg versus m H +. 



both rrih and rrin are close to the value of m#+ , which is helpful in reducing the L2HDM 
contribution to the precise electroweak data such as mw and Zfr couplings at Z-pole. For 
rriH+ > 250 GeV, the data require | sin(/3 — a)\ ~ 1 23j, and in this case, has little effects 



on the data so that it can deviate significantly from m#+. 

In summary, the above results indicate that the preferred parameter space of the L2HDM 
is 37 ^ tan (3 ^ 80, rriA 30 GeV and the other Higgs bosons lighter than 250 GeV. Note 
that the above favored region is obtained by considering all the constraints (1-11), instead 
of any individual constraint. For example, for tan (3 > 200, our results indicate that the 
CP-odd Higgs boson A as heavy as 120 GeV can still explain a M ; but such a large tan/3 is 
disfavored by the Zfr coupling at Z-pole or by r leptonic decay. Another point we should 
address is that in the L2HDM, the processes B — > V s 7 and T — > A'y cannot impose any 
further constraints The reason is that in the surviving parameter space, tan/? must 
be larger than 37 and, consequently, the couplings of bottom quark with H + and h are 
suppressed, as shown in Eq.®. Finally, we would like to emphasize that in contrast to 
the L2HDM which has a large parameter space to account for the discrepancy without 
conflicting with other experimental data, the popular type-II 2HDM is very difficult to do 
so bJ. This is one of the virtues of the L2HDM. 
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IV. IMPLICATION ON HIGGS PHENOMENOLOGY 



Eq.® indicates that the lepton couplings of A, H and H + are enhanced by large tan/3, 
while quark couplings are suppressed. Since the allowed parameter space has a large tan (3, 
the couplings of r lepton with A, H and H + are larger than the top quark couplings. So 
these scalars will decay dominantly into r leptons rather than into top quarks (if kinetically 
allowed). Moreover, a light A can change the phenomenology of other Higgs bosons by 
opening new decay modes like h, H — > AA, h,H — > AZ and H + — > AW + . As discussed 
earlier, in case of a large tan /? and a small a, /t is the SM-like Higgs boson, mainly responsible 
for the electroweak symmetry breaking and couples to weak gauge bosons like the SM Higgs. 
Therefore, the phenomenology of h is of primary importance and will be studied in the 
following. 




T 
■c 
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FIG. 3: Same as Fig.l, but showing the branching ratio of h — > AA versus m^. 



Fig. 3 shows the branching ratio of h — > AA versus m^. Here we have considered all 
the decay modes of h including h — > VV, AZ, rf, bb, ti. This figure shows that for most 
of the allowed parameter space (about 99%), h —>■ AA — > rfrf is the dominant decay 
mode. This will make the detection of h difficult at the LHC because the lightness of A 
(note rriA < 20 GeV for about 70% surviving samples) will make the r leptons from its 
decay highly collimated 0, 24], which is usually regarded as a difficult scenario in Higgs 



9 



discovery at the LHC in the next-to-minimal supersymmetric model 2j. Another reason for 
the detection difficulty of h is that for more than 80% of the allowed parameter space, the 
width of h is found to be larger than 10 GeV. Such a wide width will smear the peak of the 
invariant mass distribution of /i-decay products and make the detection more difficult. 

We note that in the L2HDM, A — > is the second largest decay mode of A. So 

h — > AA can give the multi-muon signal, like the scenario proposed in 25]. Unfortunately, 
in the L2HDM the branching ratio of A —>■ fi + fi~ is of order 10~ 3 , which will make the 
channel h — > AA — > 4/i quite hopeless at the LHC. Note that some authors have considered 



2fi + 2 jets 



261 ] as well as 



the channels h -> AA -> 2fi + 2r |22j and /i -> -> 4r 

the diffractive Higgs production pp ^ pp + h followed by h — > At [27| to detect such a h, 
but all these studies did not consider the worse case of a fat h. We also checked that the 
branching ratio of h —>■ 77 is usually suppressed to be less than 10~ 6 and thus too small for 
the detection. 




FIG. 4: Same as Fig.l, but showing the cross section of HA associated production at the LHC 
versus run- 



Furthermore, we examine other complementary new channels for detecting the Higgs 
sector of the L2HDM. Firstly, we check the associated HA production at the LHC. The 
cross section of this process is shown in Fig. 4 and one can learn that for m# < 140 GeV 
the cross section is larger than 100 fb. The dominant decay of H in this case is found to be 
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H — > AA, with a branching ratio larger than 80%, so the main signal of this process is 6r. 
Due to the lightness of A, great efforts are needed to analyze the signal and the backgrounds 
in order to detect this process at the LHC Secondly, we note that A is always lighter than 
Z boson in the allowed parameter space and thus it may be produced from Z decays. So 
we investigate the decay Z — > ttA and find its branching ratio ranging from 10~ 5 to 10~ 4 
for rriA < 40 GeV (corresponding to 98% of the allowed parameter space). Such a large rate 
is within the sensitivity of the GigaZ option at the proposed International Linear Collider 
28|. 

V. CONCLUSION 

We examined various direct and indirect constraints on the lept on- specific two-Higgs 
doublet model and then checked the property of the Higgs bosons in the allowed parameter 
space. We found that the allowed space has a very light CP-odd Higgs boson A (m^ < 30 
GeV with 95% possibility) which is composed dominantly by the leptonic Higgs and decays 
dominantly into t + t~. The SM-like Higgs boson h decays dominantly in the mode h — > 
AA — ► At, which may make the Higgs discovery difficult at the LHC. We also checked other 
possibilities for testing the Higgs sector of this model and found that the decay Z — > t + t~A 
has a branching ratio ranging from 10~ 5 to 10 -4 , which may be accessible at the the GigaZ 
option of the ILC. 
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